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1.  Introduction 

For  multi-stage  linear  programs  with  the  staircase  structure 
it  has  been  observed  that  a nested  decomposition  (Staircase)  al- 
gorithm [23, [3],  can  become  more  efficient  than  a direct  Simplex 
approach.  In  general  this  tendancy  increases  with  increasing  problem 
size.  However,  even  for  smaller  problems,  the  Staircase  algorithm 
usually  converges  rapidly  during  Phase  1 amd  the  beginning  of 
Phase  2 before  exhibiting  a long  "tail"  towards  optimality.  This 
suggests  a hybrid  algorithm  using  Staircase  initially  and  then 
switching  to  Simplex.  Since  the  Staircase  algorithm  produces 
solutions  which  are  in  general  nonbasic,  a special  interfacing 
procedure  is  required  in  Simplex  to  adopt  nonbasic  starting  solu- 
tions . 
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Computational  experience  indicates  that  the  hybrid  algorithm  j 
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*1  *1  * ‘’l 

(1) 

x^iO,  t*l,...,T 


where  c^ 

is 

1 X n^,  x^  is 

n^  X 1,  A^  is  m^  X n^,  is 

“t+i  * "t 

and 

d^  is  m^  X 1 

in  dimensions. 

In  [2]  and  [3],  the  Staircase  algorithm  is  developed  from  an 
application  of  nested  decomposition  to  (1) . The  original  problem 
is  replaced  by  a sequence  of  smaller,  independent  subproblems  co- 
ordinated through  price  amd  proposal  communication  (in  the  sense 
of  Damtzig  eind  Wolfe  [1])  between  adjacent  svibproblems . The  algo- 
rithm seeks  an  optimal  coordination  (if  one  exists)  which  in  turn 
deteraines  em  optimal  solution  to  the  original  problem. 

This  is  done  in  three  phases.  Phase  1 seeks  an  initial  feas- 
ible solution  to  (1)  or  shows  that  none  exists.  Phase  2 seeks  an 


optimal  solution  to  (1)  or  shows  that  it  is  unbounded  from  below. 

Throughout  Phase  2,  primal  feasibility  is  maintained  implicitly. 

A Phase  3 procedure  is  required  to  reconstruct  a feasible  solu- 
tion. This  is  normally  done  at  optimality.  Dual  feasibility  is 

i 

also  maintained  during  Phase  2 so  that  a lower  bovind  on  the  objec-  ; 

! 

tive  is  available  as  an  optimality  criterion. 

It  has  been  observed  in  [3]  that  relative  to  a direct  simplex 
approach  the  Staircase  algorithm  becomes  more  efficient  with 
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increasing  problem  size.  However,  the  threshold  problem  size 
differs  considerably  for  different  classes  of  Staircase  problems. 
Moreover,  even  for  smaller  problems,  the  Staircase  algorithm 
usually  converges  rapidly  during  Phase  1 and  the  beginning  of 
Phase  2 before  exhibiting  a relatively  long  "tail"  of  convergence 
towards  optimality. 

3.  A Hybrid  Algorithm 

The  edxsve  observation  suggests  using  Staircase  to  obtain  a 
near-optimal  solution  and  then  switching  over  to  Simplex.  Hence, 
a Hybrid  algorithm  will  be: 

Step  (i) : Phase  1 of  Staircase 

vintil  feasibility. 

Step  (ii) : Phase  2 of  Staircase 

until  objective  is  within  p % of 
best  lower  bound  where  p is  user  supplied. 
(p=®  implies  shipping  of  Step  (ii)). 

Step  (iii) : Phase  3 of  Staircase 

for  current  feasible  solution 

Step  (iv) ; Staircase-Simplex  Interface 

from  solution  in  Step  (iii)  to 
basic  feasible  solution. 

Step  (v) : Phase  2 of  Simplex. 

until  optimality. 
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Step  (iv)  is  necessary  because  the  solution  obtained  in 
Step  (iii)  is  in  general  nonbasic.  The  relative  amount  of  work 
in  this  extra  step  determines  how  well  the  Hybrid  algorithm  can 
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more  they  have  been  used  previously  in  ad  hoc  methods  to  partition 
structured  problems,  such  as  Staircase  problems.  The  general  idea 
is  to  manually  or  heuristically  partition  the  right  hand  side 
(resource)  vector  among  the  submodels  (time  periods  or  divisions) 
which  are  then  solved  independently.  Assuming  a feasible  partition 
was  chosen,  a non-basic  feasible  solution  will  then  be  available. 
This  non-basic  solution  (names  and  values  of  the  variables  not  at 
bound,  together  with  bound  information)  is  then  used  as  a starting 
point  for  the  entire  undecomposed  problem,  and  the  BASIC  procedure 
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is  used  to  derive  a basic  feasible  solution  with  at  least  as  good 
a value  as  the  non-basic  solution. 

Although  the  BASIC  proceedure  has  been  implemented  many  times, 
the  methodology  does  not  seem  to  have  been  published.  We  therefore 
give  an  outline  of  the  method  (there  are  several  variations) . 


Consider  the  problem  of 
minimizing  cx 

s\jb j ect  to  Ax*b 


with  X a nonbasic  feasible  initial  solution. 

Let  J = {j|Xj>0}  and  set  x ^ , jeJ  nonbasic  at  a temporary  bound  (TB) 

of  X with  appropriate  modification  of  the  right  hand  side.  Then  proceed 
as  follows: 

Step  (0)  Start  with  all  logical  (or  artificial)  basis. 

Step  (1)  Stop  if  J-j&.  Price  only  Xj,jeJ  for  reduced 

cost  d j . 

Step  (2)  If  dj  < 0,  raaTce  x^  basic  by  increasing  above 
X j , i.e.  treating  TB  as  a lower  bound.  If 
dj  2 0,  make  x^  basic  by  decreasing  below  x ^ , 
i.e.  treating  TB  as  an  upper  bound. 
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Step  (3)  Remove  TB  for  and  j from  J. 


Return  to  Step  (1) . 

Note  that  with  this  procedure  the  logical  basis  in  Step  (0) 
will  be  feasible  at  zero  value  since  x is  a feasible  solution. 


Therefore,  initially  m^ny  x^,  jeJ  will  enter  the  basis  without 


changing  their  values  from  x^.  Much  of  the  computation  in  pricing 


and  pivot  selection  can  be  avoided  if  we  first  scan  J and  pivot 


into  the  basis  any  x^  on  any  row  with  a nonzero  entry  in  the  up- 


dated Xj  column  and  a zero  entry  in  the  updated  right-hand-side. 


This  "crash"  procedure  is  equivalent  to  introducing  the  largest 


independent  subset  of  x ^ , jeJ  into  the  basis  before  setting  the 


remainder  to  temporary  bounds.  Experiments  have  shown  that  this 
is  essential  to  the  efficiency  of  the  interface. 


5.  Computational  Experience 

The  experimental  codes  are  written  in  FORTRAN  for  the  CDC 
7600  at  BrooTchaven  National  Laboratory.  SIMPLEX  is  based  on  the 
revised  simplex  code  LP Ml  with  product  form  of  inverse  (cf  [5^). 
STAIRCASE  is  the  implementation  described  in  [3^  , also  based  on 
LPMl.  HYBRID  is  based  on  STAIRCASE,  SIMPLEX  and  the  interface 
procedure  in  Section  4. 

Computational  experience  with  three  small-  to  medium-size 
problems  are  reported  in  Tables  1 and  2 and  Figure  1.  It  is  ob- 
served that  HYBRID  can  be  used  to  combine  the  advantages  of  STAIRCASE 


and  SIMPLEX. 
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6.  Conclusion 


Since  the  performance  of  special  purpose 
algorithms  is  uaually  highly  problem  dependent,  the  flex- 
ibility provided  by  a hybrid  algorithm  would  be  useful  in  a truly 
versatile  Staircase  decomposition  system.  Furthermore,  basic 
solutions  are  usually  preferable  from  a practical  point  of  view, 
since  fewer  activities  are  "active".  The  advantage  of  the  hybrid 
algorithm  we  propose  over  the  ad  hoc  use  of  partitioning  and  the 
BASIC  procedure,  is  that  the  Staircase  algorithm  appears  to  be 
successful  in  producing  good  feasible  solutions  relatively  quickly 
and  does  not  rest  on  unreliable  user  partitioning  of  the  problem. 

Our  preliminary  computational  results  would  indicate  that  the 
hybrid  approach  may  be  one  of  the  more  successful  for  the  notoriously 
stubborn  class  of  Staircase  models. 
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Figure  1.  Solution  history  of  Problem  SCFXMl. 
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